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The intruder feature of 31 Mg and 
the coexistence of many particle and many hole states 
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The low- lying level structure of 31 Mg has been investigated by the antisymmetrized molecular 
dynamics (AMD) plus generator coordinate method (GCM) with the Gogny D1S force. It is shown 
that the N=20 magic number is broken and the ground state has the pure neutron 2p3h configuration. 
The coexistence of many particle and many hole states at very low excitation energy is discussed. 

PACS numbers: Valid PACS appear here 



The breakdown of the neutron magic number N=20 
and the change of the shell order in neutron rich nu- 
clei have been one of the great interest in nuclear struc- 
ture studies. Experimentally, it was firstly pointed out 
from the observation of the anomalous binding energy 
and spin-parity of 31 Na 0, More convincing evi- 
dence was given by the observation of the small exci- 
tation energy [H, 0| and the large E2 transition prob- 
ability 0] of the 2f state of 32 Mg. Theoretically, the 
Hartree-Fock calculation [|| showed that shape transi- 
tion from spherical to prolate shapes possibly took place 
in neutron-rich N=20 nuclei. Shell model calculations 
@> H, that allowed particle hole excitations across 
the N=20 shell gap explained these abnormal properties 
and suggested strong deformation of 31 Na and neighbor- 
ing nuclei caused by the inversion between the normal 
and intruder co nfig urations. Since then, many experi- 
mental fTol [Ti l . Jl2t f]~3L Jlil [Tp . Hp . IT7I and theoretical 
[H, El, HO, HH, |22, M, M, [23 HI] studies have 
been made and now the systematic breaking of N=20 
magic number is theoretically investigated and some of 
them are observed. Thus, our knowledge of the break- 
down of the neutron magic number N=20 in the neutron 
rich nuclei has been increasing rapidly in this decade. 

To understand the shell structure and the shell order, 
the odd-neutron nuclei such as 31 Mg and 29 Ne will pro- 
vide essential information, since the property of the sin- 
gle particle orbital of the last neutron particle or hole 
will be strongly reflected to the property of total system. 
Recently the spin-parity of the ground state of 31 Mg was 
experimentally fixed to l/2 + [lfj. This result suggests 
that the ground state of 31 Mg has the intruder configu- 
ration in which two neutron occupy the pj-orbital and a 
neutron hole is in the Nilsson orbital [200,1/2]. However, 
mainly because of the limitation of the models, theoreti- 
cal investigation of odd-neutron nuclei is not as vigorous 
as that of even-even nuclei. Indeed, there are only a 
few shell model calculations QJ, [H El, El for 31 Mg, 
and they do not reproduce the ground state properties 
or they study the ground band and a few excited states 
using the modified interaction to reproduce the ground 
state properties. 

In this work, we discuss the structure of the 31 Mg based 
on the theoretical framework of AMD+GCM [2a Hi. 



Since we do not assume the time reversal symmetry and 
we perform the parity projection before the variation, we 
can appropriately investigate the positive- and negative- 
parity states and the coexistence of the many particle- 
hole states in odd-mass nuclei. It will be shown that the 
ground state of 31 Mg has the pure intruder configura- 
tion that has the 2p3h configuration with respect to the 
N=20 shell closure. We have successfully reproduced the 
observed properties of 31 Mg such as the magnetic mo- 
ment and the /3-decay strengths. We also predict that 
other many particle-hole states (Oplh, \p2h and 2>pAh) 
also coexist at very small excitation energy. 

In the AMD, the intrinsic wave function is given by a 
Slater determinant of the single particle wave packets, 
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where (pi is the ith single particle wave packet consisting 
of the spatial fa, spin Xi and isospin & parts. The local 
Gaussian located at Zj is employed as fa. 



fair) 

Xi 
6 



exp 



a=x,y,z 



■{r a 



proton or neutron. 



Here the centroids of the Gaussian Zj, the spin direc- 
tion a,i and Pi and the width parameters v x ,v y and v z 
are the variational parameters. As the variational wave 
function, we employ the parity projected wave function, 
$*■ = P^^. The Gogny D1S force is employed as 
an effective nuclear force V n . Coulomb force V c is ap- 
proximated by the sum of seven Gaussians. 

The calculation is performed in the following three 
steps. First, we carry out the energy variation and opti- 
mize the variational parameters. The energy variation is 
made under the constraint on the matter quadrupole de- 
formation [3 to obtain the optimum wave function $*' (/3) 
for given value j3. We note that in this calculation, we 
do not make any assumption on the spatial symmetry of 
the wave function and do not put any constraint on the 
quadrupole deformation parameter 7. Therefore, 7 has 
the optimal value for each given value of (3. The definition 
of (3 and 7 is given in Ref. [3l|. Then, we project out the 
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eigenstate of the total angular momentum J from the op- 



timized wave function, §{j K ((3) 



is referred as AMD wave function in the following. Fi- 
nally, we superpose all of the AMD wave functions on the 
energy surface that have the same parity and the angular 
momentum but have different K and (3, 
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The coefficients c n , c' n ,... are determined by solving 
the Hill- Wheeler equation (34[. <l>^ 7r will be referred as 
AMD+GCM wave function. 

We have investigated the single particle configuration 
of each intrinsic wave function by the AMD+HF method 
3l| . Once we obtain the intrinsic wave function, <j>; nt {(3) , 
we transform the single particle wave packet tp+ into the 
orthnormalized basis cp a , 
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where A Q and Ci a are the eigenvalue and eigenvector of 
the matrix Bij = ((pi\ipj). From this new basis function, 
we construct and diagonalize the Hartree-Fock single par- 
ticle Hamiltonian, 
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FIG. 1: Energy surfaces of 31 Mg as functions of matter 
quadrupole deformation parameter /3 for the positive- and 
negative-parity states. The estimated neutron particle hole 
configurations are also shown around the energy minima. 
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The eigenvalue and eigenvector provide the single par- 
ticle energy e s and the single particle wave function tp s 
= Y) a —i fas^fa- To investigate the properties of the sin- 
gle particle wave functions, we also calculate the mag- 
netic quantum number {ip s \j s) and the amount of the 
positive-parity component (ift s \P + \4> s ). We note that by 
this method we obtain only the occupied states. In the 
following, we denote the neutron m particle and n hole 
configuration relative to the N=20 spherical shell closure 
as mpnh. 

Fig. Q] shows the calculated energy surface of the 
positive- and negative-parity states. There are three 
local minima in the positive-parity state and two local 
minima in the negative-parity state. Though we do not 
assume the axial symmetry, the wave functions on the 
energy surfaces are almost axially symmetric except for 
those located at the transitional region between the en- 
ergy minima. These minima have different neutron single 
particle configurations. Fig. [2] shows the single particle 
orbitals occupied by the last 11 neutrons in the positive- 
parity state (a) and the negative-parity state (b). In 
this figure, the asymptotic quantum number [TV, n z , l z ,jz] 
of each orbital evaluated from its parity and magnetic 
quantum number is also shown. In the following, we 



use the asymptotic quantum number to specify each neu- 
tron orbital. Fig. [2] (a) shows that in the spherical re- 
gion of the positive-parity state, the last neutron parti- 
cle or the neutron hole occupies [2,0,2,3/2] orbital that 
originates in the 0c?3/2-orbital. Therefore this configu- 
ration is understood as the Oplh configuration. Around 
(3 ~ 0.3, the neutron configuration changes from Oplh 
to 2p3h. The last two neutron occupy [3,3,0,1/2] orbital 
that intrudes from the p/-orbital. And a neutron hole 
occupies [2,0,0,1/2] orbital that originates in the 0c? 3 / 2 - 
orbital. Around [3 ~ 0.5, the neutron configuration 
changes from 2pZh to 4p5h. [2,0,2,3/2] and [2,0,0,1/2] 
orbitals are no longer occupied. Both of [3,3,0,1/2] and 
[3,2,1,3/2] orbitals intruding from the p/-orbital are oc- 
cupied by two neutrons. A neutron hole is in [2,0,2,5/2] 
orbital that originates in the 0d5/2-orbital. In the same 
way, Fig. [5] (b) shows that there are two different neu- 
tron configurations lp2h and 3pAh in the negative-parity 
state. The lp2h configuration has the last neutron in 
[3,3,0,1/2] orbital and [2,0,2,3/2] orbital is unoccupied. 
Around [3 ~ 0.4 the configuration changes from \p2h to 
3p4/i. The 3p4h configuration has the last neutron in 
[3,2,1,3/2] orbital and both of [2,0,2,3/2] and [2,0,0,1/2] 
orbitals are unoccupied. Thus, the analysis of the single 
particle configurations has revealed that the minima of 
the positive-energy state correspond to the neutron Oplh, 
2p3h and Ap5h configurations and those of the negative- 
parity state correspond to the neutron \p2h and 3p4h 
configurations. It is also mentioned that protons are al- 
ways below the Z=20 shell closure and do not excited 
from the Op-orbital to the IsOd-orbital, though the or- 
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(a) single particle energy (positive-parity) 




(b) single particle energy (negative-parity) 
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FIG. 2: The single particle orbitals of the last 11 neutrons 
of the positive-parity state (a) and the negative-parity state 
(b). Filled (open) symbols show the orbitals occupied by two 
(one) neutrons. Circles (boxes) show the orbital in which 
the amount of the positive-parity component is larger than 
(smaller than) 50%. 



der of the proton [2,2,0,1/2] and [1,0,1,1/2] orbitals has 
changed in the strongly deformed region. Therefore, the 
parity of total system depends only on the neutron con- 
figuration. There are some additional comments on Fig. 
Q] and [U (1) All of the mpnh configurations except for 
the Apbh configuration appear at very small excitation 
energy. The energy difference between the lowest one 
(Opl/i) and the highest one (3p4h) is less than 5 MeV. 
This leads to the coexistence of these particle-hole config- 
urations within small excitation energy after the angular 
momentum projection as shown later. (2) The behavior 
of the neutron single particle orbitals of the positive- and 
negative-parity states are qualitatively similar to each 
other. As deformation becomes larger, [3,3,0,1/2] and 
[3,2,1,3/2] orbital originate in the ^/-orbital come down 
and [2,0,2,3/2] and [2,0,0,1/2] orbital originate in the 
0d 3 / 2 -orbital go up and cross to each other. The change 
of the particle hole configuration occurs just around the 
crossing points of these orbitals. This explains why the 
particle hole configurations appear in the order of the 
particle and hole numbers as deformation becomes larger. 
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FIG. 3: Level scheme of 31 Mg obtained by the AMD wave 
functions at energy minima and by the AMD+GCM calcula- 
tion together with the partial level scheme suggested in Ref. 
[Ir3 ]. Solid (dashed) lines show the states dominated by the 
neutron Oplh or lp2h (2p3h or 3p4/i) configurations. Numbers 
above or blow the states show the calculated and observed ex- 
citation energy in MeV. 



It is also mentioned that the behavior of the single par- 
ticle orbital is qualitatively understood by the Nilsson 
model and it may explain why the Nilsson model calcu- 
lation was qualitatively successful in Ref. [17[. (3) The 
energy of the Apbh configuration is exceptionally higher 
than others and it is understood as follows. As explained 
above, the change of the neutron particle-hole configu- 
ration is caused by the crossing of four neutron orbitals 
that originate in the pf- and Osld-orbitals. Therefore, 
up to four particle or four hole configuration can appear 
at small excitation energy, but more than five particle or 
hole configuration appear at rather high excitation en- 
ergy. This feature is common to other neighboring nuclei 
such as Ne isotopes (2(| |3(| and Mg isotopes [27| • 

After the variational calculation, we have performed 
the angular momentum projection of the wave functions 
on the energy surfaces. The left side of FIG. [3] shows the 
excitation energies of the AMD wave functions at the 
energy minima. In the case of the positive-parity state, 
the Oplh minimum generates the 3/2 + , 5/2 + and 7/2+ 
states and the 2p3h minimum generates the K v =l/2 + ro- 
tational band with negative decoupling parameter. The 
Apbh minimum generates the K 7T =5/2 + rotational band, 
but it is not shown since it is not involved in the low- 
lying states even after the GCM calculation. In the case 
of the negative-parity state, lp2h minimum generates the 
K v =l/2~ rotational band with the large negative decou- 
pling parameter a ~ —4 (the last neutron that occupy 
the O/7/2 gives a = —4). It is interesting that the lowest 
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state of lp2h configuration is not 7/2" but 3/2" because 
of its deformation. The 3p4h configuration generates the 
K w =3/2~ rotational band. As clearly seen, the energy 
of the deformed states becomes lower. For example, the 
3/2" state that has the 3p4h configuration is lowered by 
about 5.1 MeV from the minimum of the energy surface, 
while the 3/2+ state that has the Oplh configuration is 
lowered only by about 1.8 MeV. This leads to the co- 
existence of many particle-hole states within quite small 
excitation energy. However, it should be noted that there 
is no 1/2+ state except for the 2p3h configuration. It is 
also mentioned that the K quantum number of the de- 
formed bands (K 7r =l/2 + , 1/2" and 3/2") correspond to 
the asymptotic quantum number j z of the last neutron 
particle or hole. 

To complete our calculation, we have performed the 
GCM calculation by superposing all of the AMD wave 
functions on the energy surfaces. The calculated level 
scheme is shown in the middle of FIG. [3l There are mix- 
ing between different particle- hole configurations. For 
example, in the case of the 3/2+ states, the mixing be- 
tween the Oplh and 2p3h configurations lowers the en- 
ergy of the 3/2+ state that has the dominant 2p3>h con- 
figuration, while it pushes up the 3/2^" state that has 
the dominant Oplh configuration. Since the AMD+GCM 
wave function is the sum of many Slater determinants, 
its particle-hole configuration is ambiguous. To evaluate 
it, we have calculated the squared overlap between the 
AMD+GCM wave function of each state and the AMD 
wave function at each energy minima that has the definite 
particle-hole configuration. The energy gain by such con- 
figuration mixing does not exist in the case of the 1/2+ 
state that purely consist of 2p3h configuration (it has no 
overlap with the Oplh AMD wave function but has 93% 
overlap with the 2p3h AMD wave function), but it is the 
ground state with a small margin and consistent with the 
experiment [16j |. The calculated total binding energy is 
242.9 MeV. The excitation spectrum also shows the good 
agreement with the most recent spin-parity assignment 
1 16]. The first excited state is the 3/2+ state that has the 
dominant 2p3h configuration. The squared overlap with 
the 2p3h AMD wave function amounts to 61%. The low- 
lying negative-parity states appear as the second (3/2]"), 
third (7/2^) and fourth (3/2^") excited states with the 
dominant lp2h, lp2h and 3p4h configurations, respec- 



tively. The normal configuration (Oplh) appears as the 
fifth excited state 3/2^ at 0.81 MeV. Thus, from Oplh 
to ipAh states, many particle-hole configurations coexist 
below 1 MeV. 
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TABLE II: The electric quadrupole moment of the positive- 
parity states calculated by the AMD+GCM wave function 
(Q) and the rigid rotor approximation (Qs). 



Finally, we discuss observables and compare them with 
available experimental data. The magnetic moments 
of the ground state (1/2+) and the first excited state 
(3/2+) calculated without the spin quenching factor are 
— 0.91/i/v and 0.43/xn, respectively. The magnetic mo- 
ment of the ground state shows the reasonable agreement 
with the observed value — 0.89/xat [l6| and the shell model 
result [Uj, while 3/2+ state has the opposite sign to the 
ground state. The log ft value of the (3 decay from 31 Na 
and to 31 Al are summarized in TABLE fl] The ground 
state of 31 Na and the low-lying states of 31 A1 are also 
calculated by the AMD+GCM. The squared overlaps be- 
tween the AMD+GCM wave function and the AMD wave 
function are also listed to show the particle- hole nature of 
the low-lying states 31 Mg and 31 Al. It is found that the 
ground state of 31 Na (3/2+) is the mixture of the neutron 
OpOh, 2p2h and 4p4h configurations. The squared over- 
laps amount to 10%, 70% and 5%, respectively. These 
configurations feed the Oplh, 2pZh and Apbh configura- 
tions of 31 Mg. Here, we must recall that the ground 
state of 31 Mg purely consist of the 2p3h configuration. 
Therefore, it is populated only from the 2p2h configu- 
ration of 31 Na that is the dominant component of the 
31 Na. Indeed, the calculated log ft shows the stronger 
31 Na(g.s.) — > 31 Mg(g.s.) transition than others. It is 
also notable that the excited states of 31 Mg dominated 
by the 2p3h configuration have stronger transition than 
those dominated by the lp2h configuration. Though the 
spin-parity assignment of the excited states are not fixed 
yet, the tendency of the log ft is comparable with the 
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experiment 0, [nj. The /3-decay from 31 Al also reflects 
the specific charactor of the 31 Mg. Since, the ground 
state of 31 Mg purely consist of the neutron 2p3h con- 
figuration, it feeds only the neutron 2p4h configuration 
of 31 A1. The low-lying states of 31 A1 are dominated by 
the neutron normal configuration (0p2h) and the intruder 
configuration (2p4ft). The 1/2+, 3/2^, 3/2$, and 1/2$ 
states dominated by the neutron 2pAh configuration have 
smaller log ft than those dominated by other neutron 
configurations. The strength of log ft and the excitation 
energy shows good agreement with the experimental data 
[l7j . The spin-parity assignment of those 31 Al states are 
not experimental fixed yet, but the present result and 
the shell model calculation (TtI [25[ shows the reason- 
able agreement. Table. HI] summarizes the quadrupole 
moments of the positive-parity states calculated by the 
AMD+GCM wave function and those evaluated by the 
rigid rotor approximation [35| . 
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where the intrinsic quadrupole moment Qq is calculated 
from the intrinsic wave function $J lt (/3) that is the domi- 
nant component of each state. The K quantum number is 
taken as 1/2 for the states dominated by the \p2h config- 
uration and 3/2 for those dominated by the 3p4h config- 
uration. Due to the strong deformation, the yrast states 
that have the intruder configuration show large collectiv- 



ity. It is also notable that the rigid rotor approximation 
qualitatively reproduces these states, even though there 
are the configuration mixing in the AMD+GCM calcu- 
lation. 

To summarize, we have investigated the low-lying level 
structure of 31 Mg. The observed property of the ground 
state including the spin-parity is successfully reproduced. 
The coexistence of the spherical normal configuration 
and the strongly deformed intruder configuration at vary 
small excitation energy is predicted. The analysis of the 
single particle orbital has shown that the neutron Oplh, 
\p2h, 2p3h and 3p4h configurations appear on the en- 
ergy surface. The angular momentum projection leads 
to the coexistence of these configurations and the GCM 
calculation has shown the mixing between them. Among 
the obtained states, only the ground state consists of the 
pure neutron 2p3h configuration and its spin-parity l/2 + 
originates in a neutron hole in the Nilsson [2, 0, 0, 1/2] or- 
bital. This specific charactor is reflected to the /3-decay 
strengths from the ground state of 31 Na and to the low- 
lying states of Al. 
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